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Abstract 

A very simplified model of the Universe is considered in order to propose an al- 
ternative approach to the irreversible evolution of the Universe at very early times. 
The entropy generation at the quantum stage can be thought as a consequence of an 
instability of the system. Then particle creation arises from this instability. 
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Entropy generation in cosmological particle creation has been extensively studied dur- 
ing the last years. Essentially, the formalisms were developed in the framework of quantum 
field theory in curved spacetime and at this level it was already demonstrated the dissipa- 
tive nature of this process |l|, ^, ||, Q . The notion of entropy and its "emergence" in many 
of these descriptions [§, [| is based on the concept of "coarse graining," by means of pro- 
jection operator techniques. This approach is also successfully realized in the framework 
of quantum cosmology ||. However, entropy generation was really understood only in the 
later stages of cosmic evolution. In the quantum creation regime, from the usual point of 
view, it was not clear which mechanism leads to entropy generation. It is argued that the 
evolution of the system follows strictly quantum mechanical laws which are "time-reversal 
invariant." Then, without external influence, we cannot expect any change in entropy 

An alternative point of view is to consider that quantum mechanical laws are not 
"always" time-reversal invariant. There are many models for quantum unstable systems 
(See, e.g. Ref. @) ; among them, a very relevant case is the Friedrichs model |0], which 
was widely studied using the rigged Hilbert space formulation of quantum mechanics 
This model becomes of interest to us because our cosmological model can be reduced to 
the Friedrichs one in a short time approximation, as we shall see. 

Let us introduce the model. We shall use a quantum, conformally coupled, massive 
scalar field to represent the matter degrees of freedom. Then, the action S reads 

sTd^dy* + tm 2 + l -m§ 2 j, (i) 

where m stands for the mass of the fields quanta. 
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We shall also assume that the metric corresponds to a closed Friedmann-Robertson- 
Walker model. Hence, the scalar field must be homogeneous in order to satisfy the Einstein 
equations. Nevertheless, we want to go beyond the two degrees of freedom of this min- 
isuperspace and deal with a statistical system. So, to obtain a consistent model with 
inhomogeneous matter fields, we shall not impose the momenta constraints but only the 
Hamiltonian constraint averaged over each spatial hipersurface. Formally, it is expressed 
in the fact that the lapse function iV depends only on time (for a detailed discussion see 



where gij is the induced metric on the spatial hipersurface labeled with time t. We can 
express an arbitrary scalar field configuration in terms of the eigenfunctions of the spatial 
Laplacian, the spherical harmonics on the 3-sphere Qf m - 



Ref. §). 



Then the interval results 



ds 
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N 2 {t)dt 2 + q 2 (t)g ij dx i dx j 



(2) 



$(t,x) = -X>n(*)Qn(x) 



(3) 



n 



where n denotes the set {n,l,m}. 



Therefore, we can derive the Hamiltonian 




where ir = (—m 2 /N)q and p n = (q/N)<p n are the canonical momenta associated with q 
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and (p n , respectively. 

At this point, we shall not follow the standard canonical quantization procedure. In 
order to recover a notion of time in quantum gravity (for a valuable treatise on this 
subject see Ref. [||), we can break the temporal reparametrization invariance and choose 
a privileged time. In this framework, we can use the "probabilistic time," previously 



introduced by one of us [1C]. 



We break the temporal reparametrization invariance fixing the lapse function N. Each 
particular choice for the fixing implies a different probabilistic time. In our case, the more 
appropriated one is to choose N = q, and 9 turns out analogous to the conformal time in 



quantum field theory in curved spacetimes (for details see Refs. 12, ^3[). Then, following 



the concepts of Ref. [11], we can obtain a Schrodinger equation, where the temporal 



evolution is labeled by the probabilistic time 



= { 2 - it* 2 + \ ? + (n2 + m2q2) A k (5) 



The Hamiltonian can be rewritten in terms of the usual creation and annihilation 
operators as 



f°° + r°° 

H = -n a j a+ dnnblbn + X dn g (n)(a + o t ) 2 (6 n + %) 2 + const, (6) 
Jo " Jo 

where we have regularized the expression of H and we have slightly modified it without 
substantially affecting the qualitative behavior of the system. This trick allows us to show 
the dissipative characteristic of the system because the Poincare period becomes infinite 

mm. 
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It is very hard to deal with the interaction term. However, if we take some physical 
processes into account, we can obtain a diagonalization of the Hamiltonian. It is easy to 
verify that the interaction term does not conserve the particle number, it "connects" states 
which differ in a multiple of pairs of particles. For giving an explicit matrix representation, 
we adopt the basis associated to the "particle number representation," since in the limit of 
weak particle creation (which will be sufficient to demonstrate the existence of dissipation) 
the number of particles is an adiabatic invariant. In addition, we shall consider particle 
creation in the first stages of the cosmic evolution, when occupation numbers are small 
and spontaneous particle creation dominates over stimulated creation Then, in a 
first approximation, at small times since Universe creation, we can consider the creation 
of a pair of particles in each mode, and only in the lowest occupation number states. 
This sort of approximation is widely used in nuclear physics, where a collective state 
is obtained by diagonalizing the interaction in a limited number of shell-model states of 
particle excitation. It is usual to refer to such calculations as the "Tamm-Dancoff method" 
J15|| . Then, under these conditions we can consider only the following states: 



where the first place corresponds to the quanta associated to a, and the following places 
correspond to the bath (scalar field) modes. So, the Hamiltonian in this sector reads 



1;0,...,0,...> = |1) and |1; .., 2 n , ...) = |fi) 



3V2\g(n) 



...\ 



H = 



3V2\g(Q) 



n + 2Q + 12Xg(n) + 3XJ g(Q)dfl 



(7) 



V 



•■■/ 
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Redefining 



u>q = —Qo + 3A / g(Q)d£l, 
Jo 



(8) 
(9) 



we obtain the Hamiltonian of a Friedrichs model for this subspace (one-particle sector of 
an oscillator in a bosonic reservoir |T^1 ), in the form 



H 



^A 5 (^p) ...\ 



(10) 



It can be rewritten as 



H = U0 \l)(l\+ £ to [f + 6A, MM+A jT ^ 5 (^) (H>H+I"><1|), 

(11) 

where uiq G 7£>o and <?(0) = 0. The associated Hilbert space is H = C © £ 2 (^,oo). 

The eigenvalue problem can be exactly solved (see Ref. ||). This well-known solution 
includes eigendistributions which must be interpreted in an adequate extension of the 
Hilbert space. 

In order to make explicit an exponential decaying law we can analytically extend the 
problem to the complex plane (for a detailed derivation see Refs. |8| |12[)- In this case, 
the second-sheet extension of the reduced resolvent of the Hamiltonian, 



6 



oT l {z) 



z — loq — -A 



z _|_ 6A5( ^o; 



presents a complex pole zq. Thus, in a generalized spectral decomposition of the Hamilto- 
nian zq appears as a "complex eigenvalue" ||]. For small values of A, zq can be estimated 



MS 



UJQ+V 



duo 



|AV(^Q) 



+ i 



9 \2 2, UJ 0\ 

"V 9 ( T } 



(12) 



where V denotes principal part, zq is usually expressed as 



1 

z =uj - i- (7 > 0), 



(13) 



where u)o — is the level shift and I/7 is the mean-life of the unstable level. 

By performing complex integration and using the residues theorem, a complex contri- 
bution coming from the pole zq gives rise to a set of states which we shall use to obtain a 
spectral resolution of H. The associated eigendistributions (corresponding to the discrete 
part) from right and left, respectively are Q 



|l-) = K /+ (z )]- 1 / 2 |l)+/ r ^ 



zo-f -6Ag( - 2° ) 



(14) 



(1- 



(^o)]- 1/2 ((l|+/r^ 



fo 3 



2 -f -6Ag( — 
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where the curve T passes just under the pole zq and II stands for the second Riemann- 
sheet. The vectors corresponding to the continuous solution are 



M + 



3^2 M5£) 

2 a + (u) 



3V2 



(15) 



where a± = a(u> ± ie) and 



1 1 +27ri ^i+ 6A 5(^ a )-^o) 



The orthonormality and completeness relations read 



(1+|1-) = 1, (l + |uT) = = 0, <w + |a/ - ) = S(u-u'), 



+ l u duj\u-){uj + \ = 1. 
2 



(16) 
(17) 



Finally, the spectral decomposition of H is 



H = zo|l~)(l + | + / duu\u~)(u^ 



(18) 



The eigenvectors ( |l7| ) (from right and left, respectively) and their associated complex 
eigenvalue zq have not sense in the standard Hilbert space formulation of quantum me- 
chanics. However, H admits a self-adjoint extension to a "rigged Hilbert space," where 
the states corresponding to complex eigenvalues acquire meaning [Jsj, |l7], O]. The idea is 



to restrict the topology of the Hilbert space 7i in order to have a nuclear space ft such that 
fi C H. Then, the linear functionals on Q belong to a bigger space fi x , the topological 
dual of fi, such that we build up a Gel'fand triplet [jig ]: C H C O x . The nuclear space 
is chosen by using convergence prescriptions. In our case, this condition is fulfilled if we 



use Hardy class functions |l^, ^C§ in order to define two nuclear subspaces + and f2_, 
for the upper and the lower half-plane, respectively (Hardy classes from above and below, 
respectively) pl| . Prom the Paley- Wiener theorem [ 19 j we know that the time evolution 
splits into two semigroups, i.e. 



|l-(0)> = e- iHe \l~) = e-^° e e-i e \l-), 6 > 0, 

(19) 

|l+(0)) = e - iHe \l+) = e- i *° e el 9 \l+}, 6 < 0. 

Then, the state \1~(9)} decays towards the past and the state \l + (6)) decays towards 
the future. These vectors belong to the dual spaces of f2+ and respectively 0^ and fi* . 
The physical meaning of the f^-states is the following: the state |1~(6>)) corresponds to 
the evolution from an unstable state towards a stable one. The state |l + (#)) corresponds 
to the formation process (which is usually neglected in the case of a large meandife). 

The survival amplitude of the state 1 1} in time 9 is given by 

1 CI poo 2/ ui—WQ \ 

(lle-^ll) = , , , e~ iz ° e + -A 2 / e-^ e 9 ,\ 2 , dw. (20) 

The first term is dominant over the second for short times (with the exception of the 
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short initial Zeno period p2|| ), showing an exponential behavior. On the other hand, the 
second term is the continuous background and is relevant for later times. In this way the 
initial matter-less state |1) is unstable, favouring transitions to states with greater number 
of particles with a smooth exponential law. 

Formally, this dissipative behavior can be made explicit by the following "H-theorem." 



Let us consider the operator [23] 



y = \\-){\-\. (21) 



It evolves according to the Heisenberg equation of motion 



y{Q) = e iH0 y(O)e- iHe = e^ e y(0). (22) 



y{6) is a monotonic increasing function in time 9, i.e. 



Any component of this operator is therefore a Lyapunov function. Hence, in this model, 
"dissipation" results as a consequence of the resonance between the scale factor and the 
scalar field (this fact was already suggested as the responsible of an instability at the 
classical level with the consequent apparition of chaos |24j ) , promoting states with greater 
occupation numbers, as is usually expected. Thus, dissipation can be thought as a gen- 
eralized particle creation process, the vacuum of the matter field is unstable under pair 
creation, whereby changes in the "system" (i.e. the degree of freedom related to gravi- 
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tation, the scale factor) lead to excitations of degrees of freedom of the "bath" (i.e. the 
infinity degrees of freedom of the scalar matter field) , which then propagate or diffuse away 
Q. What we have found is that these excitations are originated in an intrinsic resonance 
(associated with the pole zq) between the system and the bath. 

Although the model we have presented here cannot be considered as a very realistic 
one, it constitutes a first step for an alternative approach to the problem of finding the 
yet "unclear" mechanism of dissipation at the early stage of the Universe. 



Acknowledgements 

This work was partially supported by Universidad de Buenos Aires. 



11 



References 

[1] B.L. Hu and D. Pavon, Phys. Lett. B180 (1986) 329. 

[2] B.L. Hu and H.E. Kandrup, Phys. Rev. D35 (1987) 1776. 

[3] H.E. Kandrup, Phys. Rev. D37 (1988) 3505 ; D47 (1988) 3413. 

[4] E. Calzetta and B.L. Hu, Phys. Rev. D40 (1989) 656. 

[5] E. Calzetta, M. Castagnino, and R. Scoccimarro, Phys. Rev. D45 (1992) 2806. 

[6] V. Weisskopf and E.P Wigner, Z. Phys. 63 (1930) 54; 65 (1930) 18; T.D. Lee, Phys. 
Rev. 95 (1954) 1329; V. Glaser and G. Kallen, Nucl. Phys. 2 (1956) 706 ; PT. 
Matthews and A. Salam, Phys. Rev. 112 (1958) 283; K. Nakanishi, Prog. Theor. 
Phys. 19 (1958) 607; M. Levy, Nuovo Cimento 13 (1959) 115; J. Schwinger, Ann. 
Phys. (NY) 9 (1960) 169; L. Fonda, G.C. Ghirardi, and A. Rimini, Rep. Prog. Phys. 
41 (1978) 587. 

[7] K.O. Friedrichs, Comm. Pure Appl. Math. 1 (1948) 361. See also E.B. Davies, J. 
Math. Phys. 15 (1974) 2036; P. Exner, Open Quantum Systems and Feynman integrals 
(Reidel, Dordrecht, 1985). 

[8] E.C.G. Sudarshan, C.B. Chiu, and V. Gorini, Phys. Rev. D18 (1978) 2914 ; G. 
Parravicini, V. Gorini, and E.C.G. Sudarshan, J. Math. Phys. 21 (1980) 2208; E.C.G. 
Sudarshan, Phys. Rev. A46 (1992) 37; E.C.G. Sudarshan and C.B. Chiu, Phys. Rev. 
D47 (1993) 2602; I.E. Antoniou and I. Prigogine, Physica A192 (1993) 443. 



12 



[9] K.V. Kuchaf , Time and interpretations of Quantum Gravity in "Proceedings of the 
4 th Canadian Conference on General Relativity and Relativistic Astrophysics," eds. 
G. Kunstatter, D. Vincent, and J. Williams, (World Scientific, Singapore, 1992). 

[10] M. Castagnino, Phys. Rev. D39 (1989) 2216. 

[11] M. Castagnino and F. Lombardo, Phys. Rev. D48 (1993) 1722. 

[12] M. Castagnino, F. Gaioli, and E. Gunzig, Cosmological features of time asymmetry, 
Preprint IAFE (1995). 

[13] M. Castagnino, F. Gaioli, and D. Sforza, Substantial irreversibility in quantum cosmol- 
ogy in "Proc. of the V Int. Workshop on Instabilities and Nonequilibrium Structures," 
Santiago de Chile, eds. E. Tirapegui and W. Zeller (Kluwer Acad. Publ., Dordrecht, 
1995). 

[14] M. Castagnino, E. Gunzig, I. Nardone, I. Prigogine, and S. Tasaki, Quantum cosmol- 
ogy and large Poincare systems in "Fundamentals papers in Theoretical Physiscs," 
ed. N. Namiki (APS Publ., New York, 1994). 

[15] A.M. Lane, Nuclear Theory: Pairing force correlations and collective motion (W.A. 
Benjamin, New York, 1964). 

[16] N. Bleistein, H. Neumann, R. Handelsman, and L.P. Horwitz, Nuovo Cimento A41 
(1977) 389. 

[17] A. B6hm, J. Math. Phys. 22 (1981) 2813; A. B6hm, M. Gadella, and G.B. Mainland, 
Am. J. Phys. 57 (1989) 1103; A. Bohm and M. Gadella, Dirac kets, Gamow vectors, 



13 



and Gel' f and triplets, "Lectures Notes in Physics," vol. 348 (Springer- Verlag, Berlin, 
1989). 

[18] I.M. Gel'fand and N.Y. Vilenkin, Generalized Functions, vol. 4 (Academic Press, New 
York, 1964); A. Bohm, The Rigged Hilbert Space and Quantum Mechanics, "Lectures 
Notes in Physics," vol. 78 (Springer- Verlag, Berlin, 1978). 

[19] See, for example, P. Duren, Theory of H p - spaces (Academic Press, New York, 1970). 

[20] M. Gadella, J. Math. Phys. 24 (1983) 1462; 24 (1983) 2142; 25 (1984) 2481. 

[21] M. Castagnino, M. Gadella, F. Gaioli, and R. Laura, Gamow vectors and time asym- 
metry, Preprint IAFE (1995). 

[22] B. Misra and E.C.G. Sudarshan, J. Math. Phys. 18 (1977) 756 ; C.B. Chiu, E.C.G. 
Sudarshan, and B. Misra, Phys. Rev. D16 (1977) 520. 

[23] H. Hasegawa, T. Petrosky, I. Prigogine, and S. Tasaki, Found. Phys. 21 (1991) 263. 

[24] E. Calzetta and C. El Hasi, Class. Quantum Grav. 10 (1993) 1825. 



14 



